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We extract parameters relevant for distinguishing among single-field inflation models from the 
Wilkinson Microwave Anisotropy Probe (WMAP) three-year data set, and also from WMAP in 
combination with the Sloan Digital Sky Survey (SDSS) galaxy power spectrum. Our analysis leads 
to the following conclusions: 1) the Harrison-Zeldovich model is consistent with both data sets at 
a 95% confidence level; 2) there is no strong evidence for running of the spectral index of scalar 
perturbations; 3) Potentials of the form V tx <f) p are consistent with the data for p = 2, and are 
marginally consistent with the WMAP data considered alone for p — 4, but ruled out by WMAP 
combined with SDSS. We perform a "Monte Carlo reconstruction" of the inflationary potential, and 
find that: 1) there is no evidence to support an observational lower bound on the amplitude of grav- 
itational waves produced during inflation; 2) models such as simple hybrid potentials which evolve 
toward an inflationary late-time attractor in the space of flow parameters are strongly disfavored by 
the data, 3) models selected with even a weak slow-roll prior strongly cluster in the region favoring 
a "red" power spectrum and no running of the spectral index, consistent with simple single-field 
inflation models. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Inflation [l| has become the dominant paradigm for un- 
derstanding the initial conditions for structure formation 
and for Cosmic Microwave Background (CMB) tempera- 
ture anisotropics. In the inflationary picture, primordial 
density and gravitational-wave fluctuations are created 
from quantum fluctuations, "redshifted" beyond the hori- 
zon during an early period of superluminal expansion of 
the universe, then "frozen" l2j.yj.ljf. Perturbations at the 
surface of last scattering are observable as temperature 
anisotropies in the CMB, as first detected by the Cosmic 
Background Explorer satellite The latest and most 

impressive confirmation of the inflationary paradigm has 
been recently provided by the three-year data from the 
Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite HHHII3]- The WMAP collaboration has produced 
new full-sky temperature maps in five frequency bands 
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from 23 to 94 GHz based on the first three years of the 
WMAP sky survey. The new maps, which are consistent 
with the first-year maps and more sensitive, incorporate 
improvements in data processing made possible by the 
additional years of data and by a more complete analy- 
sis of the polarization signal. WMAP data support the 
inflationary model as the mechanism for the generation 
of super-horizon curvature fluctuations. 

The goal of this paper is to make use of the re- 
cent WMAP three-year data (WMAP3) to discriminate 
among the various single-field inflationary models. As 
such, this paper represents a complete update of our pre- 
vious analysis of the first-year WMAP data. 

For single-field inflation models, the relevant parame- 
ter space for distinguishing among models is defined by 
the scalar spectral index n, the ratio of tensor to scalar 
fluctuations r, and the running of the scalar spectral in- 
dex dn/dhik. We employ Monte Carlo reconstruction, 
a stochastic method for "inverting" observational con- 
straints to generate an ensemble of inflationary poten- 
tials compatible with observation 0, 0] . In addition to 
encompassing a broader set of models than usually con- 
sidered (large-field, small-field, hybrid and linear mod- 
els) , Monte Carlo reconstruction makes it possible easily 
to include effects to higher order in slow roll. 
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The paper is organized as follows: In Sec. II we will 
quickly review single-field inflation models and their ob- 
servables. In Sec. Ill we define the inflationary model 
space as a function of the slow-roll parameters e and 77. 
In Sec. IV we describe our analysis method as well as our 
results. Since a study of the implications of the WMAP3 
data for single field models of inflation has been already 
performed by the WMAP collaboration themselves [3, 
we will also specify briefly the differences between our 
analysis and theirs. In Sec. V we describe a Monte Carlo 
reconstruction method to determine an ensemble of infla- 
tionary potentials compatible with observations. In Sec. 
VI we present our conclusions. 



II. SINGLE-FIELD INFLATION AND THE 
INFLATIONARY OBSERVABLES 

In this section we briefly review scalar field models of 
inflationary cosmology, and explain how we relate model 
parameters to observable quantities. Inflation, in its most 
general sense, can be defined to be a period of acceler- 
ating cosmological expansion during which the universe 
evolves toward homogeneity and flatness. Within this 
broad framework, many specific models for inflation have 
been proposed. We limit ourselves here to models with 
"normal" gravity (i.e., general relativity) and a single 
order parameter for the vacuum, described by a slowly 
rolling scalar field <fi, the inflaton. 

A scalar field in a cosmological background evolves 
with an equation of motion tfi + 3H<j) + V' (</>) = 0. The 
evolution of the scale factor is given by the scalar-field 
dominated FRW equation, 



H l = 



~<t> 2 + V(<t>) 



v(<p)-4> 2 



(i) 



We have assumed a flat Friedmann-Robertson- Walker 
metric g^ v = diag(I, — a 2 , — a 2 — a 2 ), where a 2 (t) is the 
scale factor of the universe. Inflation is defined to be a pe- 
riod of accelerated expansion, a > 0. A powerful way of 
describing the dynamics of a scalar field-dominated cos- 
mology is to express the Hubble parameter as a function 
of the field <fi, H — -ff (</>), which is consistent provided cf> 
is monotonic in time. The equations of motion become 

0EEE3: 

[H >m 2 - l ^H\4>) = -^V{4>). (2) 

These are completely equivalent to the second-order 
equation of motion. The second of the above equations 
is referred to as the Hamilton- J acobi equation, and can 
be written in the useful form 



1 - \<<t>) 



( 8tt 



V 377J 2 



PI 



V(<f>), 



(3) 



where e is defined to be 



4vr \H{(f>) 



(4) 



The physical meaning of e(<p) can be seen by expressing 
Eq. as 



H 2 (cj>) [l-e(cj>)], 



(5) 



so that the condition for inflation, (a/a) > 0, is equiva- 
lent to e < I. The scale factor is given by 



a cx e 



N 



exp 



Halt 



where the number of e-folds N is 



N 



Hdt 



H 



■d(j) 



d<j> 



mpi 7 0c y/^j' 



(0) 



(7) 



We will frequently work within the context of the slow- 
roll approximation, which is the assumption that the evo- 
lution of the field is dominated by the drag from the cos- 
mological expansion, so that <fi ~ and <j) — —V/3H. 
The equation of state of the scalar field is dominated by 
the potential, so that p ~ —p, and the expansion rate 
is approximately H 2 ~ 8ttV (cf>) / Sm 2 ^. The slow roll ap- 
proximation is consistent if both the slope and curvature 
of the potential are small, V , V" <C V. In this case the 
parameter e can be expressed in terms of the potential as 



'pi 



'pi 



(8) 



4tt V H ((f)) J 16-k V V O) 
We will also define a second "slow-roll parameter" j] by 



8tt 



H(4>) 
V(<t>) 



V(4>) 
V(4>) 



(9) 



Slow roll is then a consistent approximation for e, r\ <C I . 

Inflation models not only explain the large-scale homo- 
geneity of the universe, but also provide a mechanism for 
explaining the observed level of inhomogeneity as well. 
During inflation, quantum fluctuations on small scales 
are quickly redshifted to scales much larger than the hori- 
zon size, where they are "frozen" as perturbations in the 
background metric. The metric perturbations created 
during inflation are of two types: scalar, or curvature per- 
turbations, which couple to the stress-energy of matter 
in the universe and form the "seeds" for structure for- 
mation, and tensor, or gravitational-wave perturbations, 
which do not couple to matter. Both scalar and ten- 
sor perturbations contribute to CMB anisotropy. Scalar 
fluctuations can also be interpreted as fluctuations in the 
density of the matter in the universe. Scalar fluctuations 
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can be quantitatively characterized by the comoving cur- 
vature perturbation Pn- As long as the equation of state 
e is slowly varying, the curvature perturbation can be 
shown to be Q 



Pn /2 (*) 



(") = 


' H 1 " 






mpi y/iTe_ 


k=aH 



(10) 



The fluctuation power spectrum is in general a function 
of wavenumber k, and is evaluated when a given mode 
crosses outside the horizon during inflation, k = aH. 
Outside the horizon, modes do not evolve, so the ampli- 
tude of the mode when it crosses back inside the horizon 
during a later radiation- or matter-dominated epoch is 
just its value when it left the horizon during inflation. 
Instead of specifying the fluctuation amplitude directly 
as a function of k, it is convenient to specify it as a func- 
tion of the number of e-folds N before the end of inflation 
at which a mode crossed outside the horizon. 
The spectral index n for Pn is defined by 



1 



d\nP n 



(11) 



din A; ' 

so that a scale-invariant spectrum, in which modes have 
constant amplitude at horizon crossing, is characterized 
by n = 1. 

The power spectrum of tensor fluctuation modes is 
given by 



Pt' 2 (k N ) = 



AH 



(12) 



N 



The ratio of tensor-to-scalar modes is then Pt / Pn = 16e, 
so that tensor modes are negligible for e <C 1. 



III. THE INFLATIONARY MODEL SPACE 

To summarize the results of the previous section, infla- 
tion generates scalar (density) and tensor (gravitational- 
wave) fluctuations which are generally well approximated 
by power laws: Pn{k) oc fe n_1 , Pt (k) oc k nT . In the 
limit of slow roll, the spectral indices n and ut vary 
slowly or not at all with scale. We can write the spectral 
indices n and tit to lowest order in terms of the slow roll 
parameters e and n as: 



n • 
ut ' 



1 - 4e + 2t], 
-2e. 



(13) 



The tensor/scalar ratio is frequently expressed as a quan- 
tity r, which is conventionally normalized as 



16e 



Pn 



(14) 



1 This normalization convention is differ ent from that used in our 
analysis of the first- year WMAP data llll . which used the con- 
vention r = lOe. In this paper, we have adopted the more stan- 
dard normalization convention r = 16e. 



The tensor spectral index is not an independent param- 
eter, but is proportional to the tensor/scalar ratio, given 
to lowest order in slow roll by ut — — 2e = — r/8. This is 
known as the consistency relation for inflation. A given 
inflation model can therefore be described to lowest or- 
der in slow roll by three independent parameters, Pn, 
Pt, and n. If we wish to include higher-order effects, we 
have a fourth parameter describing the running of the 
scalar spectral index, dn/d\nk. 

Calculating the CMB fluctuations from a particular 
inflation model reduces to the following basic steps: (1) 
from the potential, calculate e and r\. (2) From e, cal- 
culate N as a function of the field (j>. (3) Invert N (</>) 
to find 4>n. (4) Calculate Pn, n, and Pp as functions of 
4>, and evaluate them at <f> — (f>N. For the remainder of 
the paper, all parameters are assumed to be evaluated at 
4> — 4>n, where N varies from 46 to 60. 

With the normalization fixed, the relevant parameter 
space for distinguishing between inflation models to low- 
est order in slow roll is then the r — n plane. (To next or- 
der in slow-roll parameters, one must introduce the run- 
ning of to.) Different classes of models are distinguished 
by the value of the second derivative of the potential, 
or, equivalently, by the relationship between the values 
of the slow-roll parameters e and 77. Each class of mod- 
els has a different relationship between r and n. For a 
more detailed discussion of these relations, the reader is 
referred to Refs. Gj|. 

Even restricting ourselves to a simple single-field infla- 
tion scenario^the number of models available to choose 
from is large [l| . It is convenient to define a general clas- 
sification scheme, or "zoology" for models of inflation. 
We divide models into three general types: large-field, 
small-field, and hybrid, with a fourth classification, lin- 
ear models, serving as a boundary between large- and 
small-field models. 

First order in e and rj is sufficiently accurate for the 
purposes of this Section, and for the remainder of this 
Section we will only work to first order. The generaliza- 
tion to higher order in slow roll will be discussed in the 
following. 



A. Large-field models: — e < n < e 

Large-field models have inflaton potentials typical of 
"chaotic" inflation scenarios [2(j , in which the scalar field 
is displaced from the minimum of the potential by an 
amount usually of order the Planck mass. Such models 
are characterized by V" ((j>) > 0, and — e < ry < e. The 
generic large-field potentials we consider are polynomial 
potentials V (<j>) = A 4 (0/^) p , and exponential potentials, 
y(0)=A 4 exp(^). 

For the case of an exponential potential, V ((f) oc 
exp(^//i), the tensor/scalar ratio r is simply related to 
the spectral index as 



r = 8(1 -to) 



(15) 
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but the slow roll parameters have no dependence on the 
number of e- folds N. 

For inflation with a polynomial potential, V (4>) oc 4> p , 
we have 



n - 1 = - 



2+p 
2N ' 



2N \p + 2 



(1-n), (16) 



so that tensor modes are large for significantly tilted spec- 
tra. 



B. Small-field models: r] < — e 

Small-field models are the type of potentials that arise 
naturally from spontaneous symmetry breaking (such as 
the original models of "new" inflation plU^ ) and from 
pseudo Nambu-Goldstone modes (natural inflation [2j|). 
The field starts from near an unstable equilibrium (taken 
to be at the origin) and rolls down the potential to a sta- 
ble minimum. Small-field models are typically character- 
ized by V" ((f)) < and 77 < — e. Typically e (and hence 
the tensor amplitude) is close to zero in small-field mod- 
els. The generic small-field potentials we consider are of 
the form V (<fr) = A 4 [1 — (<fr/ [i) p ] , which can be viewed as 
a lowest-order Taylor expansion of an arbitrary potential 
about the origin. The cases p — 2 and p > 2 have very dif- 
ferent behavior. For p — 2, n—1 — — (l/27r)(mpi//z) 2 and 
there is no dependence upon the number of e- foldings. 
On the other hand 

r = 8(1 - n)exp[-l - N (1 - n)] . (17) 

For p > 2, the scalar spectral index is 

2 fp-1 



n ~ 1 



N \p-2 



(18) 



independent of (mpi//x). Assuming /x < mp\ results in an 
upper bound on r of 



r < 



P 



N(p-2) 



8- 



Np{p-2) 



p/(p-2) 



(19) 



a second free parameter characterizing the models. Be- 
cause of this extra freedom, hybrid models fill a broad 
region in the r — n plane. For (0jv//i) ^> 1 (where 4>n 
is the value of the inflaton field when there are N e- 
foldings till the end of inflation) one recovers the same 
results as the large-field models. On the contrary, when 
(0jv/a*) 1) the dynamics are analogous to small-field 
models, except that the field is evolving toward, rather 
than away from, a dynamical fixed point. While in prin- 
ciple "hybrid" models can populate a broad region of the 
inflationary parameter space, the presence of a dynami- 
cal fixed point means that there is a simple subclass of 
hybrid models that live in a narrow band of parameter 
space along a line with r ~ 0, n > 1, and dn/dlnk ~ 0. 
We will see below that while the WMAP3 data do not 
rule out the entire region which we label here as "hybrid," 
the simplest hybrid models evolving near the dynamical 
fixed point are clearly disfavored by the data. 

An example of a model which falls into the "hybrid" 
region of the 1 — n plane away from the dynamical fixed 
point is a potential with a negative power of the scalar 
field, V (4>) = Vo [1 + a (mpi/0) p ], used in intermedi- 
ate inflation |28j| and dynamical supersymmetric inflation 
[2^| . The power spectrum is blue: the spectral index 
given by n - 1 ~ 2(p + l)/[(p + 2)(N tot - N)}, where 
iVtot is the total number of e- foldings, and the parame- 
ter r is generally negligible. However, the model exhibits 
running of the spectral index which would be potentially 
detectable by future experiments, 



dn 
dhik 



For example, for p = 2 and n = 1.2, the running is 
dn/d\nk = —0.05 [3(|. When the running is sizable, the 
tensor contribution is totally negligible, 



r«P^ /2 (n-l)( 3 P +5 )/(P+ 2 ). 



(21) 



D. Linear models: 77 



Linear models, V (</>) oc 0, live on the boundary be- 
tween large-field and small-field models, with V" (</>) = 
and r\ = — e. The spectral index and tensor/scalar ratio 
are related as 



C. Hybrid models: < e < r\ 

The hybrid scenario HE HE H3 frequently appears 
in models which incorporate supersymmetry into infla- 
tion. In a typical hybrid-inflation model, the scalar field 
responsible for inflation evolves toward a minimum with 
nonzero vacuum energy. The end of inflation arises as a 
result of instability in a second field. Such models are 
characterized by V" (0) > and < e < rj. We con- 
sider generic potentials for hybrid inflation of the form 
V((j)) = A 4 [1 + [<i>/p) p ] ■ The field value at the end of 
inflation is determined by some other physics, so there is 



r = - (1-n). 



E. Other models 



(22) 



This enumeration of models is certainly not exhaus- 
tive. There are a number of single-field models that 
do not fit well into this scheme, for example logarith- 
mic potentials V (4>) = V Q [l + (Cg 2 /87r) In (<j>/n)] typ- 
ical of supersymmetry [jj, where C counts the degrees 
of freedom coupled to the inflaton field and g is a cou- 
pling constant. For this kind of potential, one gets 
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n-1 ~ -(l/A^)(l + 3C 5 2 /167r 2 ) andr ~ (Cg 2 /n 2 )(l/N). 
This model requires an auxiliary field to end inflation and 
is more properly categorized as a hybrid model, but falls 
into the small-field region of the r — n plane. 



F. Beyond first order 

The four classes of inflation models, categorized by the 
relationship between the slow-roll parameters as — e < 
1) < e (large field), r\ < — e (small field, linear), and 
< e < 77 (hybrid), cover the entire r — n plane and 
are in that sense complete at first order in the slow roll 
parameters. However, this feature is lost going beyond 
first order: models can evolve from one region to another. 
This feature is manifest when changing the parameter TV, 
and is particularly relevant for those models for which 
the running of the observables with the scale is sizable 
[3l| . Therefore, it is important to realize that the lowest- 
order correspondence between the slow-roll parameters 
and the class of models does not always survive to higher 
order in slow roll. For instance, for potentials of the 
form V ((f)) = A 4 / ((f) I p) , the parameter A is generally 
fixed by CMB normalization, leaving the mass scale /z 
and the number of e-folds N as free parameters. For 
some choices of potential, for example V oc exp ((f) /n) or 

V oc 1 — (4>/fi) 2 , the spectral index n varies as a func- 
tion of /1. These models therefore appear for fixed N as 
lines on r — n plane. Changing N results in a broadening 
of the lines. For other choices of potential, for example 

V oc f — ((f) I ' n) p with p > 2, the spectral index is indepen- 
dent of /1, and each choice of p describes a point on the 
zoo plot for fixed N. A change in N turns each of these 
points into lines, which smear together into a continuous 
region. 



IV. ANALYSIS AND RESULTS 

The method we adopt is based on the publicly available 
Markov Chain Monte Carlo (MCMC) package cosmomc 
[3^ . We sample the following eight-dimensional set of 
cosmological parameters, adopting flat priors on them: 
the physical baryon and CDM densities, Ub = flbh 2 and 
lu c = fl c h 2 , the ratio of the sound horizon to the angular 
diameter distance at decoupling, 9 S , the scalar spectral 
index, its running and the overall normalization of the 
spectrum, n, dn/dhxk and A at k = 0.002 Mpc -1 , the 
tensor contribution r, and, finally, the optical depth to 
reionization, r. Furthermore, we consider purely adia- 
batic initial conditions, we impose flatness, and we use 
the inflation consistency relation to fix the value of the 
tensor spectral index tit- We also restrict our analysis to 
the case of three massless neutrino families; introducing 
a neutrino mass in agreement with current neutrino os- 
cillation data doesn't change our results in a significant 
way. 

We include the three-year data Q (temperature and 



polarization) with the routine for computing the likeli- 
hood supplied by the WMAP team and available at the 
LAMBDA web site. 2 We marginalize over the amplitude of 
the Sunyaev-Zel'dovich signal, but the effect is small: in- 
cluding/excluding the correction changes our conclusions 
on the best fit value of any single parameter by less than 
2%, and always well within the 68% C.L. contours. We 
treat beam errors with the highest possible accuracy (see 
Ref. 0, Appendix A. 2), using full off-diagonal temper- 
ature covariance matrix, Gaussian plus lognormal likeli- 
hood, and fixed fiducial C/s. The MCMC convergence 
diagnostic is done on 8 chains though the Gelman and 
Rubin "variance of chain mean" /"mean of chain vari- 
ances" R statistic for each parameter. Our f — D and 
2 — D constraints are obtained after marginalization over 
the remaining "nuisance" parameters, again using the 
programs included in the cosmomc package. In addition 
to the CMB data, we also consider the constraints on 
the real-space power spectrum of galaxies from the Sloan 
Digital Sky Survey (SDSS) 0. We restrict the analy- 
sis to a range of scales over which the fluctuations are 
assumed to be in the linear regime (k < 0.2/i~ 1 Mpc). 
When combining the matter power spectrum with CMB 
data, we marginalize over a bias b considered as an 
additional nuisance parameter. Furthermore, we make 
use of the HST measurement of the Hubble parameter 
H = lOO/i km s^Mpc -1 by multiplying the likeli- 
hood by a Gaussian likelihood function centered around 
h = 0.72 and with a standard deviation a = 0.08. Fi- 
nally, we include a top-hat prior on the age of the uni- 
verse: 10 < t Q < 20 Gyrs. 



A. Results 

As now common practice, we plot the likelihood con- 
tours obtained from our analysis on three different planes, 
r vs. n, dn/d\nk vs. n, and r vs. dn/dhik: we do so in 
Figs. dim Presenting our results on these planes is useful 
for understanding the effects of theoretical assumptions 
and/or external priors. 

We consider two different choices of datasets: the 
WMAP3 dataset alone, and WMAP3 plus the additional 
information of the SDSS. By analyzing these different 
datasets we can check the consistency of the SDSS large- 
scale structure data with WMAP3, something that is not 
completely trivial since the WMAP3 data seems to prefer 
models with a lower value for the ag, parameter than the 
one inferred from the SDSS data (see Refs. ;9] and [33jp. 

In Fig. ^ we show the 68% and 95% likelihood con- 
tours on the r vs. n plane in the case of WMAP 3 
only (left panel) and WMAP3+SDSS (right panel). We 
also consider a prior on the running: the results on 
the top panel are obtained allowing the possibility of 



2 http://lambda.gsfc.nasa.gov/ 
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FIG. 1: Constraints on the n — r plane for different choices of 
experimental datasets. The analyses in the top panels include 
a running spectral index, while the analyses in the bottom 
panels are without running. The shaded regions indicate 68% 
and 95% C.L. 



dn/dlnk ^ while on the bottom panel assume no run- 
ning. The WMAP3 only case including running exhibited 
relatively poor convergence due to a degeneracy in the 
four-dimensional parameter space of r, n, dn/dhik, and 
normalization. Adding the SDSS data set removed the 
degeneracy and substantially improved the convergence 
of the MCMC code. 

Let's first investigate the case of no running. Marginal- 
izing over all the remaining nuisance parameters we con- 
strain n and r to 0.94 < n < 1.04 and r < 0.60 at 
95% C.L. Models with n < 0.9 are therefore ruled out 
at high significance, as are models with n > 1.05. The 
data clearly set interesting constraints on tensor modes. 
Models with n < 1 must have r < 0.4 at 95% C.L. Mod- 
els with n < 0.9 must have a negligible tensor compo- 
nent. Including the SDSS data further reduces the limit 
on the amplitude of the gravitational wave component 
with a relatively smaller effect on the spectral index pa- 
rameter. For WMAP3+SDSS we constrain n and r to 
0.93 < n < 1.01 and r < 0.31. 

If we allow running the main effect is to open the con- 
tours toward higher value of n and r. With running, 
marginalizing over all the remaining nuisance parame- 
ters, we constrain n and r to 1.02 < n < 1.38 and r < 
1.09 at 95% C.L. for WMAP3 alone and 0.97 < n < 1.21 
and r < 0.38 in the case of WMAP3 plus SDSS. 

Models with n = 1 are therefore in very good agree- 
ment with CMB data in the presence of a tensor compo- 
nent or running different from zero. Of particular inter- 
est is the Harrison-Zel'dovich (HZ) model: n = 1, r = 0, 



TABLE I: One-dimensional confidence limits on inflation- 
ary parameters, marginalized over all other parameters, for 
WMAP3 alone and WMAP3 + SDSS. 



no running/ 
running 


limits on n, r 
95% C.L. 


data 
set 




0.94 < n < 1.04 
r < 0.60 


WMAP3 ONLY 
WMAP3 ONLY 


no running 


0.93 < n < 1.01 
r < 0.31 


WMAP3 + SDSS 
WMAP3 + SDSS 




1.02 < n < 1.38 
r < 1.09 
-0.17 < dn/dlnk < -0.02 


WMAP3 ONLY 
WMAP3 ONLY 
WMAP3 ONLY 


running 


0.97 <n< 1.21 
r < 0.38 
-0.13 < dn/dlnk < 0.007 


WMAP3 + SDSS 
WMAP3 + SDSS 
WMAP3 + SDSS 



dn/dlnk = 0. As we see from the bottom panel of Fig. 
n pure HZ spectra are not ruled out at more than 95% 
C.L. from CMB data alone. In particular, we found that, 
considering the whole sets of models in our 8-D chain, the 
HZ best-fit model is at Ax 2 / 2 = 2 - 04 ; 2 - 77 > and 3 -96 with 
respect to the the best fit in the case of no running and no 
tensors, including tensors but no running and including 
tensors and running. When we include the SDSS data 
we found that the HZ best fit model is at Ax 2 /2 = 3.07 
with respect to the best fit in the case of no running and 
no tensor, A% 2 /2 = 3.4 with respect to the best fit with 
no running and A% 2 /2 = 5.1 with respect to the overall 
best fit. Since Ax 2 /2 = 6.4 at 95.4% confidence level for 
6 degrees of freedom, those numbers clearly indicate that 
even when the SDSS data is included, the HZ spectrum 
is in reasonable agreement with the data. 

The fact that the scale-invariant value n = 1 is consis- 
tent with the data at the 95% C.L. when no running is 
imposed, considerably weakens the bounds on inflation- 
ary models found in Ref. [35| where the original WMAP3 
error bars were adopted concluding that n = 1 was ruled 
out at more than 99% C.L. 

In Fig.[21and Fig. 0a degeneracy is evident: an increase 
in the spectral index n is equivalent to a negative scale de- 
pendence (dn/dlrtk < 0). We emphasize, however, that 
this behavior depends strictly on the position of the pivot 
scale fc(j: choosing fco = 0.05/i Mpc -1 would change the 
direction of the degeneracy. Models with n ~ 1.1 need 
a negative running at more than about the 95% C.L. 
(about 4cr in the case of WMAP3+SDSS). It is interest- 
ing also to note that models with a red spectral index, 
n < 1.0, are in better agreement with the data with a 
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FIG. 2: Constraints on the n — dn/d\nk plane for different 
choices of experimental datasets. The shaded regions indicate 
68% and 95% C.L. 
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FIG. 3: Constraints on the dn/dlnk — r plane for different 
choices of experimental datasets. The shaded regions indicate 
68% and 95% C.L. 




1.05 



FIG. 4: The n,r parameter space WMAP3 alone (open con- 
tours) and WMAP3 + SDSS (filled contours), with a prior 
of dn/dlnk = 0. The line segments show the predictions 
for V{<j>) = m 2 4> 2 and V{4>) = X(j> 4 for N = [46,60]. The 
dashed lines show the 68% C.L. and 95% C.L. contours from 
the chains made public by the WMAP team, which do not 
include an HST prior on Ho or an age prior. The scale of 
the plot is chosen to allow direct comparison with Fig. 14 of 
Spergel et al. jg. The shaded regions indicate 68% and 95% 
C.L. 



current data, therefore, do not suggest the presence of 
running at more than 95% C.L. 

Finally, we compare our results with those presented in 
Spergel, et al. Q. While there is qualitatively good agree- 
ment, a tension appears when we compare our contour 
plots in Fig. 0] (the no running case) with those presented 
in Fig. 14 of Ref . J\ . Models with a pure HZ spectrum ap- 
pear to be ruled out by WMAP3 alone at about the 99% 
C.L. in Ref. Q, while our analysis indicates broader con- 
tours, with the Harrison-Zel'dovich spectrum inside the 
95% C.L. region. Similarly, the contours in Ref. ap- 
pear to rule out V(4>) = A0 4 , while our analysis indicates 
that this potential is still marginally consistent with the 
WMAP3 data alone at 95% confidence. In order to bet- 
ter understand this discrepancy, we compared our results 
directly with the chains made public by the WMAP team 
and available at the Lambda web site. 3 We found that the 
error contours derived from the publicly available chains 
arc considerably larger than those shown in Fig. 14 of 
Ref. 0: error contours from our analysis of the WMAP 
chains arc plotted as dashed lines in Fig. 0] 4 None of 
the contours are as tight as those shown in Spergel et al., 
and the discrepancy is significant enough to influence im- 
portant conclusions about the model space, in particular, 



zero or positive running (see Fig. l|2J), while models with 
a sizable gravity wave background need a negative run- 
ning (see Fig. EJ. For the WMAP3 alone case the run- 
ning is bounded by -0.02 > dn/d\nk > -0.17 at 95% 
C.L. (0.007 > dn/d\nk > -0.13 for WMAP3+SDSS). 
We found that the best fit from WMAP3 alone with 
dn/d\nk = is at Ax 2 / 2 = 1.2 (Ax 2 / 2 = 0.2 when 
including SDSS) with respect to the overall best fit. The 



l! http://lambda.gsfc.nasa.gov. 

4 The difference between the WMAP-team contours and our con- 
tours as plotted in Fig. [3] can be accounted for by the fact that, 
unlike the WMAP-team analysis, we include priors on Ho from 
the HST Key Project data and a top-hat age prior. We have in- 
dependently reproduced the dashed-line contours shown in Fig. 
Hi with our own analysis. 
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the consistency of a Harrison-Zel'dovich spectrum with 
the data. There appears to be a clear inconsistency be- 
tween our results and contours shown in Spergel et at, 
Figs. 12 and 14. 



V. MONTE CARLO RECONSTRUCTION 

In this section we describe Monte Carlo reconstruc- 
tion, a stochastic method for "inverting" observational 
constraints to determine an ensemble of inflationary po- 
tentials compatible with observation. The method is de- 
scribed in more detail in Refs. [l2l . In addition to 
encompassing a broader set of models than we consid- 
ered in the previous section, Monte Carlo reconstruction 
allows us easily to incorporate constraints on the running 
of the spectral index dn/d In k as well as to include effects 
to higher order in slow roll. 

We have defined the slow-roll parameters e and n in 
terms of the Hubble parameter H ((f)) in a previous sec- 
tion. Taking higher derivatives of H with respect to the 
field, we can define an infinite hierarchy of slow roll pa- 
rameters 361: 



Ah 



TOpi 
7T 

m ' 



pi 

47T 



{H'f- 1 S e+V < 



H 



(23) 



1 



Here we have chosen the parameter a = 1r\ — 4e ~ n 
to make comparison with observation convenient. 

It is convenient to use N as the measure of time during 
inflation. As above, we take t e and <p e to be the time and 
field value at end of inflation. Therefore, N is defined 
as the number of e-folds before the end of inflation, and 
increases as one goes backward in time (dt > => dN < 
0): 



d 

dN 



din a 2y / 7r 



mpi _ d 



(24) 



where we have chosen the sign convention that yft has 
the same sign as H' : 



>fe = + 



7W PL H' 

H ' 



(25) 



Then e itself can be expressed in terms of H and N simply 

as 



f dH 

HdN 



(26) 



Similarly, the evolution of the higher-order parameters 
during inflation is determined by a set of "flow" equations 



de 

In 



e (<r + 2e) . 



da 

dN 
e A H ) 



dN 



-5ecr - 12e 2 + 2 ( 2 A H ) , 
t-1 



-cr+ [i- 2)e 



A H + £+i A H .(27) 



The derivative of a slow roll parameter at a given order 
is higher order in slow roll. 

A boundary condition can be specified at any point in 
the inflationary evolution by selecting a set of parame- 
ters e, a, 2 Ah, • • • for a given value of N. This is sufficient 
to specify a "path" in the inflationary parameter space 
that specifies the background evolution of the spacetime. 
Taken to infinite order, this set of equations completely 
specifies the cosmological evolution, up to the normal- 
ization of the Hubble parameter H. Furthermore, such 
a specification is exact, with no assumption of slow roll 
necessary. In practice, we must truncate the expansion at 
finite order by assuming that the 1 Ah are all zero above 
some fixed value of I. We choose initial values for the 
parameters at random from the following ranges: 



N 
e 
a 
2 A H 
3 A H 



M+1 A 



[46, 60] 

[0,0.8] 

[-0.5,0.5] 

[-0.05,0.05] 

[-0.025,0.025], 

0. 



(28) 



Here the expansion is truncated to order M by setting 
m+1 Ah = 0. In this case, we still generate an exact solu- 
tion of the background equations, albeit one chosen from 
a subset of the complete space of models. This is equiv- 
alent to placing constraints on the form of the potential 
V ((p) , but the constraints can be made arbitrarily weak 
by evaluating the expansion to higher order. For the pur- 
poses of this analysis, we choose M = 5. The results are 
not sensitive to either the choice of order M (as long as 
it is large enough) or to the specific ranges from which 
the initial parameters are chosen. 

Solutions to the truncated flow equations are solutions 
for which all of the derivatives of the Hubble constant 
above order M + 1 vanish: 



d e H 



= 0, I > M + 2, 



with a simple polynomial solution |39| . 
H((f>)=H (l+Ai</>+--- + A M +ic 



(29) 



(30) 



The Hamilton- Jacobi Equation © can be applied to 
this solution to derive an analytic form for the poten- 
tial in terms of the parameters A\ , . . . , Am+i ■ The set of 
boundary conditions in Eq. I)28|l then consist of a weak 
slow- roll prior on the polynomial fit for H(<f>): the in- 
flaton must be slowly rolling at least at one point in its 
evolution. Thus, while the flow equations in and of them- 
selves simply define an expansion in H(<fi), the choice of 
boundary condition and the requirement that inflation 
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last at least 46 e-folds comprise a well-defined physical 
prior on the inflationary model space. 

Some interesting recent papers have explored alterna- 
tive methods for constraining the "model space" of in- 
flation. In particular, Ref. |40| incorporates the lowest- 
order flow parameters directly into the Monte Carlo 
Markov Chain fit, although they do not include effects to 
higher order in slow roll. Refs. [4lL apply a Bayesian 
model selection approach to the problem, but also do not 
consider higher-order effects which in principle contribute 
to a running spectral index. Our analysis extends these 
results by including running of the spectral index as well 
as effects to higher order in slow roll. 

Once we obtain a solution to the flow equations 
[e(N),a(N), e Xu(N)], we can calculate the predicted val- 
ues of the tensor/scalar ratio r, the spectral index n, and 
the "running" of the spectral index dn/dlnk. To low- 
est order, the relationship between the slow roll param- 
eters and the observables is especially simple: r = I6e, 
n — 1 = a, and dn/dlnk = 0. To second order in slow 
roll, the observables are given by [36L |43| , 

r = 16e [1 - C (a + 2e)] (31) 

for the tensor/scalar ratio, and 

n - 1 = a - (5 - 3C) e 2 - J (3 - 5C) ae+ - (3 - C) ( 2 A H ) 

(32) 

for the spectral index. The constant C = 4(ln2+7) — 5 = 
0.0814514, where 7 ~ 0.577 is Euler's constant. Deriva- 
tives with respect to wavenumber k can be expressed in 
terms of derivatives with respect to N as |44J 



— = - fl - ) d 
dN ( e> din A;' 



(33) 



The scale dependence of n is then given by the simple 
expression 



dn 
din k 



1 



dn 



1-eJ dN' 



(34) 



which can be evaluated by using Eq. 1-iLl and the flow 
equations. For example, for the case of V oc 4 , the 
observables to lowest order are 

16 



n- 1 

dn 
dhik 



N + l 
3 



N + r 
3 



N(N + 1) 



(35) 



The final result following the evaluation of a particular 
path in the M-dimensional "slow-roll space" is a point in 
"observable parameter space," i.e., (r,n,dn/dlnk), cor- 
responding to the observational prediction for that par- 
ticular model. 

The reconstruction method works as follows: 



1. Specify a "window" of parameter space: e.g., cen- 
tral values for n— 1, r, or dn/dlnk and their asso- 
ciated error bars. 

2. Select a random point in slow roll space, [e, n, 1 Xb], 
truncated at order M in the slow roll expansion. 

3. Evolve forward in time (dN < 0) until either (a) 
inflation ends (e > 1), or (b) the evolution reaches 
a late-time fixed point (e = Ah = 0, a = const.). 

4. If the evolution reaches a late-time fixed point, cal- 
culate the observables r, n — 1 , and dn/dhik at this 
point. 

5. If inflation ends, evaluate the flow equations back- 
ward N e-folds from the end of inflation. Calculate 
the observable parameters at that point. 

6. If the observable parameters lie within the speci- 
fied window of parameter space, compute the po- 
tential and add this model to the ensemble of "re- 
constructed" potentials. 

7. Repeat steps 2 through 6 until the desired number 
of models have been found. 

We performed the Monte Carlo reconstruction using 
the more restrictive of the data sets considered, combin- 
ing the WMAP3 data with the Sloan Digital Sky Survey. 
We ran the reconstruction code long enough (10,703,502 
iterations) to collect 10,000 models consistent with the 
WMAP3 + SDSS error bars: 4115 are within the 68% 
C.L. contours, and 5885 are within the 95% C.L. con- 
tours. 

To illustrate the degree of overlap between the various 
classes of model, the predictions for different models are 
shown in the top panel of Fig. [5] including the effect of 
the uncertainty in the number of e-folds TV. The different 
classes of potential do not have significant overlap, and 
it is therefore possible to distinguish one from another 
observationally. 

Figure also shows the points generated by Monte 
Carlo Reconstruction in the n, r parameter space. Since 
there is no measure on the space of initial conditions, the 
distribution of points generated by the flow equations 
cannot be interpreted in a rigorously statistical fashion: 
the error bars are those generated from the data using 
COSMOMC, and the points plotted are points gener- 
ated by the flow equations consistent with those errors, 
including running of the spectral index as a parameter. 
The clustering of the models in the parameter space, 
however, is significant: selecting models based on even 
an extremely weak assumption of slow-roll results in a 
strong clustering of the models in the region favoring a 
red spectrum and dn/dlnk = 0. 

In this sense, the preference for running and a blue 
spectrum present in the data itself contains very little 
information relevant to constraining slow-roll inflation 
models: it can be interpreted simply an artifact of an 
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FIG. 5: Inflationary models plotted against the 68% and 95% 
WMAP3 + SDSS error contours. Top panel: the predictions 
of various specific inflationary potentials (solid bands) plot- 
ted against the error bars from WMAP3 + SDSS with a prior 
of dn/dlnk — 0. Bottom panel: 10,000 models generated 
by flow Monte Carlo consistent with the WMAP3 + SDSS 
data sets including running as a parameter, indicated by the 
larger error contours. The contours with a dn/dhik — 
prior are plotted as a reference, and were not used in the 
Monte Carlo reconstruction. (Some data points fall outside 
the error contours plotted because likelihoods for the mod- 
els were calculated using the full three-dimensional likelihood 
function C(n, r, dn/d In k) , and the contours were obtained by 
marginalizing over dn/d In k). 



"accidental" parameter degeneracy in the data. Allow- 
ing running as a parameter but assuming slow roll in- 
flation gives constraints on the inflationary model space 
largely consistent with an analysis which assumes neg- 
ligible running as a prior on the parameter space from 
the beginning. In other words: there is no evidence for 
inflation with a measurable running of the spectral index. 

From the flow equations 1)27(1 it is evident that the line 
along the r = axis, with e = £ Ah = is a fixed point of 
the flow evolution, including taking the flow equations to 
infinite order. 5 For parameters on the "red" side of scale 
invariance, i.e. a < 0, this fixed point is unstable: flow 
moves away from the fixed point as N — > 0, and toward 
the fixed point as N — > oo. Conversely, the fixed point 
for a > is stable: models evolve toward this fixed point 
at late times, N — ► 0. Integrating the flow equations for- 
ward in time yields one of two possible outcomes. One 



See Refs. |l2tl45j for a detailed discussion of the fixed-point struc- 
ture of the slow roll space. 



possibility is that the condition e = 1 may be satisfied for 
some finite value of N, which defines the end of inflation. 
We identify this point as N = so that the primordial 
fluctuations are actually generated when N = [46,60]. 
Alternatively, the solution can evolve toward an infla- 
tionary fixed point with r — and n > 1, in which case 
inflation never stops. In reality, inflation must stop at 
some point, presumably via some sort of instability, such 
as the "hybrid" inflation mechanism H3, [2^, EE Ea • Ex- 
amples of potentials which fall into this class of models 
are the simplest hybrid potentials, 



V(4>)=A 4 







1 + 









(36) 



Here we take the observables for such models to be 
the values at the late-time attractor. Since models on 
the attractor are by definition those for which the vari- 
ation in the slow roll parameters with N vanishes, such 
models also predict zero running of the scalar spectral 
index. We find that the WMAP3 data strongly disfa- 
vor models which evolve to a late-time asymptote with 
r = 0, n > 1, and dn/dlnk = 0. The 95% confidence 
limit for a blue spectrum with no tensors and no running 
(i.e., not marginalized over r) from WMAP3 alone is 
n < 1.0007, and from WMAP3 + SDSS is n < 1.001. Of 
more than ten million models tested, only one model con- 
sistent with the data relaxed to the late-time asymptote, 
with a spectral index n = 1.0004 and r = 0.0000002; for 
all intents and purposes a Harrison-Zel'dovich spectrum. 
Every other model in the Monte Carlo reconstruction set 
was of the "nontrivial" type, with inflation ending natu- 
rally by evolving through e = 1 at late times. We note 
that the level of running required to accommodate a blue 
spectrum is severe: even dynamical supersymmetric in- 
flation, which predicts a blue spectrum and negative run- 
ning [Eq. I|20(l ]. does not produce a strong enough run- 
ning to match the data, and is also ruled out to more 
than 95% C.L. by WMAP3 + SDSS for n > 1.001. 

We can also place constraints on the energy scales rel- 
evant to inflation, in particular the "height" of the po- 
tential V ~ A 4 , and the "width" of the potential, typi- 
cally quantified as the field variation A(f> during inflation. 
Given a path in the slow roll parameter sp ace, the form of 
the potential is fixed, up to normalization |l3ll46ll47Ll48| . 
The starting point is the Hamilton-Jacobi equation, 



V(<f>) 



/3m 2 



pi 



V 87 



H 2 (<f>) 



1-^(0) 



(37) 



We have e(N) trivially from the flow equations. In order 
to calculate the potential, we need to determine H(N) 
and (f>(N). With e known, H(N) can be determined by 
inverting the definition of e, Eq. I|26|) . Similarly, (j>(N) 
follows from the first Hamilton-Jacobi equation, Eq. J2Jl: 



m PL 



dN 2,/tF 



(38) 



Using these equations and Eq. I|37|l . the form of the po- 
tential can then be fully reconstructed from the numerical 
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FIG. 6: Potentials generated by Monte Carlo Reconstruction 
consistent with WMAP3 + SDSS to 68% C.L. in the light 
shaded (yellow) region and 95% C.L. in the darker shaded 
(cyan) region. The WMAP3 data place an upper limit of 
about 2 x 10 16 GeV on the energy scale of inflation. No lower 
limit is possible without a detection of a tensor mode sig- 
nal. The concave-up line at the bottom of the figure is the 
single model in 10 7 models generated which converged to an 
inflationary fixed point at late time. 



solution for e(N) . The only necessary observational input 
is the normalization of the Hubble parameter H, which 
enters the above equations as an integration constant. 
Here we use the simple condition that the density fluctu- 
ation amplitude (as determined by a first-order slow roll 
expression) be of order 10~ 5 , 



5 P ~ H 1 
p m P i a/tte 



10" 



(39) 



A more sophisticated treatment would perform a full nor- 
malization to the CMB data |49l IHcj . The value of the 
field, (j), also contains an arbitrary additive constant. 
Fig.[S]shows the reconstructed potentials consistent with 
the WMAP3 + SDSS data set. 

We see that the energy scale for inflation favored by 
the flow Monte Carlo gives Vq between 5 x 10 14 GeV 
and 2 x 10 16 GeV, although without a detection of a 
nonzero tensor/scalar ratio, it is not possible to put a 
purely observational lower limit on the height of the in- 
flationary potential. Inflationary potentials with very 
low energy scales (in particular those with A<j) <C TOpi 
during inflation) require the imposition of a s y mmetr y 
to suppress the mass term for the inflaton [5l], 1521 l53| . 
Such potentials (corresponding to the region labeled 
V{4>) oc 1 — (0//i) p in Fig. |3J) are consistent with the 
WMAP3 data and predict an unobservably small value 
for the tensor/scalar ratio r. Even if one wished to con- 
sider such models fine-tuned (see, e.g. Ref. [M]) and 
therefore disfavored, the range of energy scales favored 
by the flow Monte Carlo (Fig. is consistent with ten- 
sor/scalar ratios as low as r ~ 10 -5 , a level unlikely to be 
detectable by any currently foreseen experiments. Refs. 
[55L |5(| suggest that fine-tuning considerations force the 



tensor/scalar ratio for a red spectrum to detectable levels 
of r ~ 0.01. We see no evidence of such an effect in the 
flow analysis, for which no explicit tuning of the potential 
is performed. 



VI. CONCLUSIONS 

In this paper, we presented an analysis of the recent 
WMAP three-year data set with an emphasis on param- 
eters relevant for distinguishing among the various pos- 
sible models for inflation. Our results are in g ood agree- 
ment with other analyses of the data [57L l58| . but show 
significant inconsistencies with the results reported by 
the WMAP team in Figs. 12 and 14 of Ref. Q. 

We found that the WMAP3 data alone are consis- 
tent within 95% C.L. with a scale-invariant power spec- 
trum, n = 1, with no running of the spectral index, 
dn/dhik — and no tensor component. The Harrison- 
Zel'dovich spectrum is therefore still not ruled out at high 
significance, a conclusion in accord with Refs. [ill l59|. 
While a detection of a running spectral index would be of 
great significance for inflationary model building [6tll6]| . 
no clear evidence for running is present in the WMAP3 
data. The data are, however, consistent with strongly 
negative running combined with a large tensor/scalar ra- 
tio and a "blue" power spectrum. 

The inclusion of the Sloan Digital Sky Survey datasets 
in the analysis has the effect of reducing the error 
bars and gives a better determination of the inflation- 
ary parameters. For instance, the inclusion of SDSS 
rules out quartic chaotic models of inflation of the form 
V((j)) ~ A0 4 . Chaotic inflation with a quadratic potential 
V(<p) ~ m 2 (f> 2 is consistent with all datasets considered. 

In addition, we applied the Monte Carlo reconstruc- 
tion technique to generate an ensemble of inflationary 
potentials consistent with observation. Our results may 
be summarized as follows: Models which evolve to a 
late-time fixed point in the space of flow parameters are 
strongly disfavored by the data. Of more than 10 million 
models analyzed in the flow Monte Carlo, one evolved to 
a late-time inflationary asymptote indistinguishable from 
a Harrison-Zel'dovich spectrum. The rest were character- 
ized by a dynamical end to inflation, with the first slow- 
roll parameter e evolving to unity in finite time. The 
late-time attractor in flow space corresponds to models 
with a "blue" power spectrum (n > 1) and negligible 
r and dn/d\nk, and we conclude that such models are 
inconsistent with the data for n > 1.001. This is a sig- 
nificant constraint on the inflationary model space. In 
particular, the data rule out the simplest models of hy- 
brid inflation of the form V{<f>) = Vo + m 2 (f) 2 as well as 
models such as V{4>) cx 1 + {fJ-/<f) F , which predict some 
negative running. This of course does not rule out all 
models for which inflation ends via a hybrid mechanism. 
Some hybrid models are characterized by a red spectrum, 
for example "inverted" hybrid models and models with 
logarithmic potentials inspired by (global) supersymme- 
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try. Finally, we find that there is no evidence to support 
any lower bound on the amplitude of gravitational waves. 
Tensor/scalar ratios as low as r ~ 10~ 5 were produced 
by the flow Monte Carlo without explicit tuning of the 
inflationary potential. 
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